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Abstract. We establish global existence and uniqueness theorems for the 
two-dimensional non-diffusive Boussincsq system with viscosity only in the 
horizontal direction, which arises in Ocean dynamics. This work improves the 
global wcll-posedncss results established recently by R. Danchin and M. Paicu 
for the Boussinesq system with anisotropic viscosity and zero diffusion. Al- 
though we follow some of their ideas, in proving the uniqueness result, we have 
used an alternative approach by writing the transported temperature (density) 
as 9 = A£ and adapting the techniques of V. Yudovich for the 2D incompress- 
ible Eulcr equations. This new idea allows us to establish uniqueness results 
with fewer assumptions on the initial data for the transported quantity 8. Fur- 
thermore, this new technique allows us to establish uniqueness results without 
having to resort to the paraproduct calculus of J. Bony. 

We also propose an inviscid o-regularization for the two-dimensional invis- 
cid, non-diffusive Boussinesq system of equations, which we call the Boussinesq- 
Voigt equations. Global regularity of this system is established. Moreover, we 
establish the convergence of solutions of the Boussinesq- Voigt model to the 
corresponding solutions of the two-dimensional Boussincsq system of equations 
for inviscid flow without heat (density) diffusion on the interval of existence 
of the latter. Furthermore, we derive a criterion for finite-time blow-up of the 
solutions to the inviscid, non-diffusive 2D Boussinesq system based on this 
inviscid Voigt rcgularization. Finally, we propose a Voigt-a rcgularization for 
the inviscid 3D Boussincsq equations with diffusion, and prove its global well- 
posedness. ft is worth mentioning that our results are also valid in the presence 
of the /3-plane approximation of the Coriolis force. 



1. Introduction 

The <i-dimensional Boussinesq system of ocean and atmosphere dynamics (with- 
out rotation) in a domain ft C K d over the time interval [0, T] is given by 

d 

(1.1a) dtu + ^djiu^u) = -Vp + 6>e d + z/Au, infix[0,T], 

(1.1b) V-u = 0, inOx[0,T], 

(1.1c) d t 6 + V ■ (u0) = kA6, inOx[0,T], 

(l.ld) u(x,0) =u (x), 0(x,O) = o (x), infi, 
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with appropriate boundary conditions (discussed below). Here v > is the fluid 
viscosity, k > is the diffusion coefficient. The spatial variable is denoted x = 
(x 1 , . . . , x d ) £ ft, and the unknowns are the fluid velocity field u = u(x, t) = 
(m 1 (x, t), . . . , u d (x, t)), the fluid pressure p(x, t), and the function 8 = 9(x, t), which 
may be interpreted physically as a thermal variable (e.g., when k > 0), or a density 
variable (e.g., when k = 0). We write e<j = (0, . . . , 0, 1) for the d th standard basis 
vector in M. d . We use the notation P° K , for the Boussincsq system with viscosity 
v > and with diffusion n > 0. We attach a subscript x to the viscosity v when we 
mean that the viscosity occurs in the horizontal direction only, i.e. in the case of 
anisotropic viscosity (see equation (1.3) below). The superscript of zero is reserved 
for a parameter a, introduced below. 

In two dimensions, the global regularity in time of the problem P® K is well- 
known (see, e.g., [8, 51]), and follows essentially from the classical methods for 
Navier-Stokes equations. However, in the case v — 0, k = 0, (Pqq), global existence 
and uniqueness still remains an open problem (see, e.g., [14, 15] for studies in this 
direction). The local existence and uniqueness of classical solutions to Pg was 
established in [15], assuming the initial data (u ,8 ) € H 3 x H 3 . In particular, an 
analogous Bealc-Kato-Majda criterion for blow-up of smooth solutions is established 
in [15] for the inviscid, non-diffusive Boussincsq system; namely, that the smooth 
solution exists on [0,T] if and only if J Q T ||V0(f)||£°o dt < oo. 

One of our main results in this study, discussed in Section 4, involves the global 
existence and uniqueness theorems for the two-dimensional non-diffusive Boussi- 
nesq system with viscosity only in the horizontal direction, denoted as P° (see 
equations (1.3) below). These equations are sometimes called the non-diffusive 
Boussincsq equations with anisotropic viscosity. In order to set the main ideas 
of our proof, in Section 3 we first establish the global existence of a certain class 
of weak solutions and the global existence and uniqueness to the two-dimensional 
viscous and non-diffusive Boussinesq system of equations (denoted as P® ) with 
Yudovich-type initial data. The other main result we have in this study is presented 
in Section 5. We propose an inviscid a-regularization for the two-dimensional invis- 
cid, non-diffusive Boussinesq system of equations (denoted as Pq ), which we call 
the Boussinesq- Voigt equations, and also establish its global regularity. We include 
in this section a study of the behavior of solutions to Pq as the parameter a — > 0, 
which leads to a new criterion for the finite-time blow-up of solutions to the 2D, or 
3D, inviscid, non-diffusive Boussinesq equations. We also give a short discussion of 
a Voigt-regularization for the three-dimensional Boussincsq equations in the case 



The two-dimensional viscous, non-diffusive Boussinesq system, (P° ) is given 




by: 
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(1.2a) 





(1.2b) 
(1.2c) 
(1.2d) 



Vu 



0. 

o, 

u (x), 0(x,O) = 0o(x) 



in T 2 x [0,T], 



d t 9 + V • (u0) 
u(x,0) 
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It has been shown in [14, 26] that the system P® , in the case of whole space 
M 2 , admits a unique global solution provided the initial data (uo,0o) G H m (M. 2 ) x 
iP™(R 2 ) with to > 3, to an integer. In fact, in [26], the authors only required 
(u , O ) G P m (R 2 ) x H m - 1 (R 2 ) with m > 3. In [14], it is shown that a Beale-Kato- 
Majda-type criterion is satisfied for the partially viscous system and therefore the 
system is globally well-posed. In [14], it is shown that the problems Pq k and P° 
both admit a unique global solution provided the initial data (uo, 0o) G P ln (R 2 ) x 
P m (R 2 ) with to > 3. Similar results are shown in [26] for Pq K but with initial data 
(u o ,0 o ) G P m (R 2 ) x i7 m - 1 (K 2 ). Global well-posedness results for rough initial 
data (in Besov spaces) is established in [25]. 

We establish in Section 3 the global well-posedness of P° in a periodic domain 
T 2 = [0, l] 2 assuming weaker initial data, namely, Uo G i? 1 (T 2 ), (we always assume 
V • u = 0) and 9 G L 2 (T 2 ). Our key idea in proving the uniqueness result is by 
writing 9 = A£, with J J2 £dx = 0, for some £, and then adapting the techniques 
of Yudovich in [27] (see also [39]). We note that the authors in [18] have shown 
the global well-posedness results in the whole space under a weaker assumption 
that Uo,6>o G L 2 (M. 2 ). The proof of their main results arise under the Besov and 
Lorentz space setting and involves the use of Littlewood-Paley decomposition and 
paradiffcrcntial calculus introduced by J. Bony [6]. We include in this study global 
well-posedness results for the problem P° under a stronger assumption on the 
initial data, namely Uo G P X (T 2 ), and 9q G L 2 (T 2 ) but using only elementary 
techniques in PDEs. Although this particular result is not an improvement to that 
of [18], we will see that applying our method in the case of anisotropic viscosity, 
we can establish an improvement to the global well-posedness results established in 
[17]. 

In Section 4, we then consider the case where the viscosity v occurs in the horizon- 
tal direction only. More precisely, assuming initial vorticity uj a G \L (defined below 
in (1.4)), initial temperature (density) 9 G L°°(T 2 ), and J J2 dx = J J2 9 dx = 0, 

we establish global well-posedness for the following system, which we denote as 
pO . 

2 

(1.3a) d t u + ^2dj(u j u) = vd\w- Vp + 9e 2 , inT 2 x[0,P], 

(1.3b) V-u = 0, in T 2 x [0,P], 

(1.3c) d t 9 + V • (u0) = 0, in T 2 x [0,r], 

(1.3d) u(x,0) = u (x), 0(x,O) = o (x), in T 2 . 

Recently, in [17], a global well-posedness result for the system P° ( m the 
whole space K 2 ), under various regularity conditions on initial data, was successfully 
established. More precisely, given that 9 G H S (R 2 ) n L°°(IR 2 ), with s G (1/2, 1], 
u G P 1 (M 2 ) and uj G L P (R 2 ) for all 2 < p < oo, and such the ujq satisfy 

(1.4) IMI^^sup L' <oo, 

p>2 VP - t 

the Boussinesq system (1.3) in the whole space with anisotropic viscosity admits a 
unique global regular solution. The condition 9o G H s with s G (5, 1] was needed for 
establishing uniqueness in [17]. We relax this condition in our current contribution. 
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We remark again that the main idea is to write 6 = A£, and then proceed using 
the techniques of Yudovich [27] for the 2D incompressible Euler equations to prove 
uniqueness. Furthermore, our method uses more elementary tools than those used 
in [17]. It is worth mentioning that very recently, in [2], the global regularity of 
classical solutions to the two-dimensional Boussinesq system in the case of vertical 
viscosity and vertical thermal diffusion was established provided an additional extra 
thermal fractional diffusion of the form (— A) d for 5 > is added. 
Let us denote by P" K the following system: 



(1.5a) -a 2 AS t u + 5 t u + ^a i (u J 'u) = -Vp + 9e d + vAu, in T d x [0,T], 

(1.5b) V-u = 0, in T d x [0,T], 

(1.5c) d t 6 + V ■ (u0) = kA6 inT d x[0,T], 

(1.5d) u(x,0) = u (x), 0(x,O) = o (x), in T l . 

In Section 5, we study in dimension d = 2 the inviscid (y — 0), Voigt-a (with 
a > 0) regularized momentum equation, namely the system Pqq, and in dimension 
d = 3 the system Pq k (with k > 0). In the case d — 2 we establish global well- 
posedness results for the problem Pq given initial data Uo G H 2 (T 2 ) with V • Uo, 
and 6q G L 2 (T 2 ). This result also hold in the easier cases k > or v > 0. In the case 
d = 3, we require n > to establish global well-poscdness results. We show that the 
problem Pq k with given initial data Uo G H 3 (T 3 ) with V • Uo, and 9$ G L°°(T 3 ) is 
well-posed globally in time. This result also hold in the easier case v > 0. Observe 
that the system Pq formally coincides with the inviscid, non-diffusive Boussinesq 
equations when a = 0. This type of inviscid a-regularization can be traced back 
to the work of Cao, et. al. [9] who proposed the inviscid simplified Bardina model 
(studied in [35]) as regularization of the 3D Euler equations. The model consists of 
the Euler equations with the term — a 2 A9 t u added to the momentum equation. We 
refer to this term as the Voigt term, and we refer to equations with this additional 
term as Voigt-rcgularized equations. The reason for this terminology is that if 
one adds the Voigt term to the Navier-Stokcs equations, the resulting equations 
happen to coincide with equations governing certain visco-elastic fluids known as 
Kelvin- Voigt fluids, which were first introduced and studied in the context of the 
3D Navier-Stokes equations by A. P. Oskolkov [40, 41], and were studied later in 
[28]. These equations are known as the Navier-Stokes- Voigt equations. They were 
first proposed in [9] as a regularization for either the Navier-Stokes (for v > 0) or 
Euler (for v = 0) equations, for small values of the regularization parameter a. 

We briefly discuss the merits of the Navier-Stokes- Voigt equations, as they are 
a special case of the Boussinesq- Voigt equations. Voigt-regularizations of parabolic 
equations are a special case of pseudoparabolic equations, that is, equations of the 
form Mut + Nu = /, where M and N are (possibly non-linear, or even non-local) 
operators. For more about pseudoparabolic equations, see, e.g., [5, 11, 19, 42, 45- 
49]. Whether in the presence of either periodic boundary conditions or physical 
boundary conditions (under the assumption of the no-slip boundary conditions 
u\qci = 0), the Navier-Stokes- Voigt equations enjoy global well-poscdness, even in 
three-dimensions), as it has been pointed out in [9]. The Euler- Voigt equations 
enjoy global well-posedness in the case of periodic boundary conditions (see, e.g., 
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[9, 34]). It is worth mentioning that the long-term dynamics and estimates for the 
global attractor, and the Gevrey regularity of solutions on the global attractor, of 
the three-dimensional Navier-Stokes-Voigt model were studied in [30] and [29], re- 
spectively. Moreover, it was shown recently in [43] that the statistical solutions (i.e., 
invariant probability measures) of the three-dimensional Navier-Stokes-Voigt equa- 
tions converge, in a suitable sense, to a corresponding statistical solution (invariant 
probability measure) of the three-dimensional Navier-Stokes equations. 

In the context of numerical computations, the Navier-Stokes-Voigt system ap- 
pears to have less stiffness than the Navier-Stokes system (see, e.g., [21, 36]). In 
[36], the statistical properties of the Navier-Stokes-Voigt model were investigated 
numerically in the context of the Sabra shell phenomenological model of turbulence 
and were compared with the corresponding Navier-Stokes shell model. 

Due to its simplicity, the Voigt a-regularization is also well-suited to being ap- 
plied to other hydrodynamic models, such as the two-dimensional surface quasi- 
gcostrophic equations, demonstrated in [31], and the three-dimensional magncto- 
hydrodynamic (MHD) equations, demonstrated in [34]. See also [21] for the appli- 
cation of Navier-Stokes-Voigt model in image inpainting. It is also worth mention- 
ing that in the case of the inviscid Burgers equation, ut + uu x = 0, this type of 
regularization leads to — a 2 u xx t + u t + uu x = 0, which is the well-known Bcnjamin- 
Bona-Mahony equation of water waves [4] . One goal of the present work is to lay 
some of the mathematical groundwork necessary to extend the Voigt regulariza- 
tion to the two-dimensional Boussinesq-equations, for the purpose of simplifying 
numerical simulations of the solutions to these equations. 

It is worth mentioning that all the results reported here are equally valid in the 
presence of the Coriolis rotation term. 



In this section, we introduce some preliminary material and notations which are 
commonly used in the mathematical study of fluids, in particular in the study of 
the Navier-Stokes equations (NSE). For a more detailed discussion of these topics, 
we refer to [16, 23, 50, 52]. 

Let T be the set of all trigonometric polynomials with periodic domain T d := 
[0, l] d . We define the space of smooth functions which incorporates the divergence- 
free and zero-average condition to be 



For the majority of this work, we take d = 2. 

We denote by L p 7 W s *, H s = W s < 2 , C°>~< the usual Lcbcsgue, Sobolcv, and 
Holder spaces, and define H and V to be the closures of V in L 2 and H 1 respectively. 
We restrict ourselves to finding solutions whose average over the periodic box T d 
is zero. Observe from the evolution equation of 9 in the Boussinesq system of 
equations (as well as the Boussinesq- Voigt system of equations), if we assume that 
the average J Td 6o(x)dx = 0, then the average of J Jd 9(x,t) dx = for all t > 0, 
and alsoj^ u(x,t) dx = for all t > provided J Td u (x) dx = 0. Therefore, we 
can work in the spaces defined above consistently. The notation V s := H s (T d ) n V 
will be convenient. When necessary, we write the components of a vector y as y J , 
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j = 1,2. We define the inner products on H and V respectively by 
(u,v) = / 



/ uV dx and ((u, v))= / djtfdjV 1 dx, 

and the associated norms |u| = (u, u) 1 / 2 , ||u|| = ((u, u)) 1 / 2 . (We use these nota- 
tions indiscriminately for both scalars and vectors, which should not be a source 
of confusion). Note that ((■, ■)) is a norm due to the Poincare inequality, (2.16), 
below. We denote by V the dual space of V. The action of V on V is denoted by 
(•, •) = (•, -) v ,. Note that we have the continuous embeddings 

(2.1) V^H^V'. 

Moreover, by the ReUich-Kondrachov Compactness Theorem (see, e.g., [1, 22]), 
these embeddings are compact. 

Following [17], we define the spaces 

: = Mini vr < °°} > 

where || • \\^ is dchncd by (1.4). This space arises naturally, due to the following 
inequality, proven in [38] (see also [17]), which is valid in two dimensions: 



(2.2) ||w|| p <CvV^lH| fflj 

for all w G i/ 1 (T 2 ), for any p G [2, oo), and where we denote by || • || p the usual L p 
norm. Note that clearly L°° C \[L C LP for every p € [2, oo). We also recall the 
well-known elliptic estimate, due to the Biot-Savart law for an incompressible vector 
field u, satisfying V ■ u = 0, and V x u = u, by means of the Caldcron-Zygmund 
theory for singular integrals: 

(2.3) !|Vuj| p <C P |M| p 

for any p G (1, oo) (sec, e.g., [27]). 

Let Y be a Banach space. We denote by L p ([0, T],Y) (which we also denote 
as LjY x ), the space of (Bochner) measurable functions t H> w(t), where w(t) G Y 
for a.e. t G [0, T], such that the integral J Q \\w(t)\\y dt is finite (see, e.g., [1]). A 
similar convention is used in the notation C k ([0,T],X) for fc-times differentiable 
functions of time on the interval [0,T] with values in Y. Abusing notation slightly, 
we write w(-) for the map t <— > w(t). In the same vein, we often write the vector- 
valued function w(-,t) as w(t) when w is a function of x and t. Wc denote by 
C°°(T 2 x [0,T]) the set of infinitely differentiable functions in the variable x and t 
which arc periodic in x with J T2 <p(-,t) dx = 0. Similarly, we denote by L P (T 2 ) = 
{<peLP(T 2 ) : J J2 ip{x) dx = 0}. 

We denote by P a : L 2 — > H the Leray-Helmholtz projection operator and define 
the Stokes operator A := -P a & with domain V(A) := H 2 n V. For <p G T>(A), 
we have the norm equivalence \Aip\ = ||<^||h-2 (see, e.g., [16, 52]). In particular, the 
Stokes operator A can be extended as a linear operator from V into V associated 
with the bilinear form ((u,v)), 

(Au, v) = ((u, v)) for all v G V. 

It is known that A -1 : H — > T>(A) >• H is a positive-definite, self-adjoint, compact 
operator from H into itself, and therefore it has an orthonormal basis of positive 
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eigenvectors {^k}^ = i in H corresponding to a non-increasing sequence of eigen- 
values (see, e.g., [16, 50]). The vectors {wsj^^ are also the eigenvectors of A. 
Since the corresponding eigenvalues of A -1 can be ordered in a decreasing order, 
we can label the eigenvalues A^ of A so that < Ai < X2 < A3 < ■ • • . Let 
H n := spanjwi, . . . , w„}, and let P n : H —> H n be the L 2 orthogonal projection 
onto H n . Notice that in the case of periodic boundary conditions in the torus T 2 
we have Ai = (2ir)~ 2 . We will abuse notation slightly and also use P„ in the scalar 
case for the corresponding projection onto eigenfunctions of —A, but this should 
not be a source of confusion. Furthermore, in our case it is known that A = — A 
due to the periodic boundary conditions (see, e.g., [16, 50]) and the eigenvectors 
Wj are of the form ake 27rlk x , with aj. • k = 0. 

It will be convenient to use the standard notation of the Navicr-Stokcs bilinear 
term 

d 

(2.4) S(wi,w a ) —P^dj 

3=1 

for Wi,W2 G V. We list some important properties of B which can be found for 
example in [16, 23, 50, 52]. 

Lemma 2.1. The operator B defined in (2.4) is a bilinear form which can be 
extended as a continuous map B : V x V — > V 1 such that 



(2.5) 



(B(wi,w 2 ),w 3 ) = / (wi • Vw 2 ) • w 3 dx, 



for every Wi,W2, W3 € V. satisfying the following properties: 
(i) For wi, W2, W3 £ V, 

(2.6) (B(w 1 ,w 2 ),w 3 ) v , = - (5(w 1 ,w 3 ),w 2 ) v , , 
and therefore 

(2.7) (B(wi,w 2 ),w 2 ) v , =0. 
(ii) For Wi, W2, W3 G V. 

(2.8) I (S(w 1 ,w 2 ),w 3 ) v , I < C|w 1 | 1 / 2 ||w 1 || 1 / 2 ||w 2 |||w 3 r /2 ||w3|| 1/2 

(2.9) I (S(w 1 ,w 2 ),w 3 ) v , I < C|w 1 | 1 / 2 ||w 1 || 1 / 2 |w 2 | 1 / 2 ||w 2 || 1 / 2 ||w3||. 

Let us define another very similar bilinear operator motivated by the transport 
term in the temperature equation. 

d 

(2.10) B(w,yj):=J2dj(w^) 

3=1 

for w G V and ip G J- with J Td ip dx = 0. We have the following similar properties 
for B which can be proven easily as in the proof of Lemma 2.1. 

Lemma 2.2. The operator B defined in (2.10) is a bilinear form which can be 
extended as a continuous map B : V x H 1 — > H , such that 



(2.11) (B(w,tp),(t)} H -i = - w-\?<t>ipdx, 

for every w G V and <fi, ip G C 1 . Moreover, 

(2.12) (B(w, ^))A)h-- = - ( fi ( w ' <t>l *)h-i . 



8 



ADAM LARIOS, EVELYN LUNASIN, AND EDRISS S. TITI 



and therefore 

(2.13) (B(w,<f>),<t>) H ^ =0. 

Furthermore, B is also a bilinear form which can be extended as a continuous map 
B : V{A) xL 2 4 H~ x . 

Here and below, C, Cj, etc. denote generic constants which may change from 
line to line. C a ,C(- • •), etc. denote generic constants which depend only upon 
the indicated parameters. K,Kj, etc. denote constants which depend on norms 
of initial data, and also may vary from line to line. Next, we recall that for an 
integrable function / such that J T2 f dx = 0, we have in two dimensions, 

(2-14) < |/| 1/2 H/|| 1/2 . 

We also recall Agmon's inequality in two dimensions (see, e.g., [3, 16]). For w G 
T> (A) we have 

(2.15) ||w|| £ - < C\^\ 1/2 \Aw\^ 2 . 
Furthermore, for all <p G V, we have the Poincare inequality 

(2.16) |M| L2 < A7 1/2 ||V^|| L2 . 

We will also make use of the following inequality, valid in two dimensions, which is 
based on the Brezis-Gallouet inequality, and which we prove in the appendix. For 
every e > 0, sufficiently small, and w G H 2 (T 2 ), 

(2.17) \\w\\ L oo < C ( IMIe- 1 / 4 + lAwIe- 1 /' 174 ' 



where C is independent of e. Finally, we note a result of deRham [52, 53], which 
states that if g is a locally integrable function (or more generally, a distribution), 
we have 

(2.18) g = Vp for some distribution p iff (g, w) = for all w G V, 

which one uses to recover the pressure. 

3. Global Well-posedness Results for the Viscous and Non-diffusive 

Boussinesq Equations. (P° ) 

Let us first define the weak formulation of problem P° K in T 2 x [0, T]. By 
choosing a suitable phase space which incorporates the divergence free condition 
of the Boussinesq equations, we can eliminate the pressure from the equation, as 
is standard in the theory of the Navier-Stokes equations. Consider the scalar test 
functions <p(x,t) G c7°°(T 2 x [0,T]), such that <p(x,T) = 0; and the vector test 
functions $(x,t) G [C 00 ^ 2 x [0,T])] 2 such that V • $(•,*) = and $(x,T) = 0. 
Then the weak formulation of problem in T 2 x [0, T] (and similarly of problem 
P° , when n = 0, in T 2 x [0,T]) is written as follows: 
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T pT 2 pT 

(u(s),&(s))ds + v ((u(s),$(s)))ds + V / (uju,0j$)ds 

Jo j=l"' 

(3.1a) = (uofa), 0)) + / (0(s)e 2 , *(*)) ds 

Jo 

(6(s),<p'(s))ds+ [ (u6, V<p) ds + k [ ((0(s), ¥>(«))) ds 
Jo Jo 

(3.1b) =(0 O (*),¥>M))- 

Remark 3.1. Note that it will become clear later that (3.1) will hold for a larger 
class of test functions, and consequently it will be sufficient to consider only test 
functions of the form 

(3.2a) $(x,t) = r m (t)e 2 " m x , with T m G [C°°([0, T])] 2 and m • T m (t) = 0, 
and 

(3.2b) V (x, t) = Xm (t)e 2mm - X , with Xm e C°°([0, T]), 

for m G (Z\{0}) 2 , since such functions form a basis for the corresponding larger 
spaces of test functions. 

In the two-dimensional case, the global well-posedness of system P° K in (1.1), 
that is, in the case k > 0, v > 0, is well-known, and can be proved in a similar 
manner following the work of [24] (see also [8, 51]). We have the following existence 
and uniqueness results for the system P° K , which will be used to prove the existence 
of weak solutions for the system P® . From here on, we only work on spaces of 
functions which are periodic and with spatial average zero. Therefore, to simplify 
notation, we write L 2 as L 2 , C as C '", etc. 

Theorem 3.2. Let T > 0, v > be fixed but arbitrary. Then, the following results 
hold: 

(i) If Uo G H, 9q G L 2 then for each k > 0, (1.1) has a unique solution 
(u K ,6 K ) in the sense of (3.1) such that u K G C([0, T],H)f] L 2 ([0, T], V), 6 K G 
C w ([0, T], L 2 ). Furthermore, there exists a constant Kq > independent of k 
such that the following bounds hold: ||u k ||£2([ ,t],v) < Kq, \\vlk\\l=°([q,t},h) < 

Kq, ||^U re || L 2 ( [ 0)T ] )V ,) < K 0) ||0re||i«>([o,T],I,*) < \6 \, ||^^ k ||l2([0,T],H- 2 ) < K 
and VK\\8 K \\L2([o,T],m) < K . 
(ii) If the initial data Uo G V and 9q G L , then the solution u K G C([0,T], V) n 
L 2 ([0, T], T>(A)) and we also have the following bounds: Iju^H i 2 ([o,T , ].X5( J 4)) < 

Ko, \\*Ik\\l°°([0,T],V) < Ko, \\foUn\\L*([0,r\,H) < K Q and \\ fa \\l*([0,T\,H-1-) < 

K . 

(Hi) 7/6*0 € L°° andu G H, then \\6 K \\L°°([o t T\,L°°) < ||0o||oo< 

(iv) If u Q G H 3 and 9 G H 2 then for each K > 0, (1.1) has a unique solution 
u K G C([0,T],# 3 ) ni 2 ([0, T],iJ 4 ) and9 K G C([0, T],H 2 ) n L 2 ([0, T],H 3 ). 

Proof. Parts (i) and (ii) are essentially proven in [8, 24, 51] following the classi- 
cal theory of Navier-Stokes equations. The uniform bounds in part (ii) will be 
established explicitly in the later proofs when called for. Part (iii) can be proven 





10 



ADAM LARIOS, EVELYN LUNASIN, AND EDRISS S. TITI 



using maximum principle and is proven for example in [8, 51]. An explicit proof 
of this theorem will also be provided below. Part (iv) can be proved using basic 
energy estimates and Gronwall's inequality again following the classical theory of 
the Navier- Stokes equations. □ 

For the current study, we now define what we mean by weak solutions and strong 
solutions for the viscous non-diffusive Boussinesq equations (P. *? ). We then state 
and prove our main results. 

Definition 3.3 (Weak solution). Let T > 0. Suppose u € H and 9 a e L 2 . We say 

that (u, 9) is a weak solution to P® Q (that is, (1.2) with k = 0) on the interval [0, T], 
if (u, 9) satisfies the weak formulation (3.1) (with k = 0), and u € L 2 ([0,T], V) n 
C([0,T],H), ^ G L^&TIV), with 9 e C([0,T],L 2 ) and f e L 1 ([0, T],H~ 2 ). 

Definition 3.4 (Strong solution). Let T > 0. Suppose #o € £ 2 and Uo € V. 
We say that (u,#) is a strong solution to P° (that is, (1.2) with k = 0) on the 
interval [0, T], if it is a weak solution in the sense of Definition 3.3, and furthermore, 
ue L 2 ([0,T], 25(A)) nC([0,T],n f € ^([0, T], ff), and ^ S L 1 ([0,T], H -1 ). 

We now state and prove our main results regarding global existence of weak and 
strong solutions to problem P° . 

Theorem 3.5 (Existence of weak solutions). Let T > be given. Let Uo e H 
and 9q € L 2 . Then there exists a weak solution of (1.1) on the interval [0,T]. 
Furthermore, system (3.1) with k = is equivalent to the functional form 

(3.3a) + vAu + B(u,u) = P a (9e 2 ) in L 2 ([0, T],V') and 

rff) 

(3.3b) — +6(u,0)=O l n i 2 ([0,T],i/- 2 ). 

at 

Moreover, if we assume 9 E L°° , then 9 e L°°([0, T],L°°). 

Proof. Our method of proof involves passing to the limit of the weak solution of 

(1.1) as k — > 0, that is, we consider n > to be a regularization parameter to 
system (1.2). Without loss of generality, we can assume < k < 1. In accordance 
with Remark 3.1 and Definition 3.3, we only consider test functions of the form 

(3.2) . 

We will show that the weak formulation 

T (u K ( S ), F^( S ) e 2 ^ m ' x ) ds + v f T ((u K (s), F m ( S )e 2 " m ' x )) ds 
o Jo 

2 r T 

f V / «u K ,r m ( s )a,e 2 " mx )d s 

3=1 Jo 

(u , F m (0)e 27 " m - x ) + / (0 K (s)e 2 , r m ( S )e 27rim ' x ) ds 
Jo 

T (9 K (s), e 2mm -*) X ' m (s) ds + [ T (u K (s)8 K (s),Ve 2 ™ m -* Xm (s)) ds 
o Jo 

(3.4b) +k f ((e K (s), e 2 ^ m ' x y m ( S ))) ds = (6 , e 2mm *) Xm (0) 
Jo 



(3.4a) 
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converges to the weak formulation of P° (see (3.1) with k = 0) as k —> 0. After 
passing to the limit in the system we then show that the limiting functions satisfy 
the aforementioned regularity properties. We proceed with the following steps. 
Step 1: Using compactness arguments to prove convergence of a subsequence. 

From Theorem 3.2, in particular from the uniform bounds (with respect to k) of 
u K , ^f-, K and ^f- in the corresponding norms, one can use the Banach-Alaoglu 
Theorem and the Aubin Compactness theorem (see, e.g., [16, Lemma 8.2] or [52]) 
to justify that one can extract a subsequence of (u K ,6 K ) (which we still write as 
(u re , K )) as k — V and elements u and 9, such that 



(3.5a) 


u« 


u 


strongly in L 2 ([0,T],H), 




(3.5b) 


u„ 


— 1 u 


weakly in L 2 ([0,T], V) and weak-* 


in L°°([0,T],H), 


(3.5c) 


du K 

dt 


du 

" ~dt 


weakly in L 2 ([0, T],V'), 




(3.5d) 


e K 


e 

_^ d9 
~dl 


weakly in L 2 ([0,T],L 2 ) and weak-* 


in L°°([0,T],L 2 ) 


(3.5e) 


dt 


weakly in L 2 ([0, T],H~ 2 ). 





Step 2: Passing to the limit in the system. 

The results from Step 1 imply that for the linear terms in (3.4), we have, by the 
weak convergence in (3.5b) and (3.5d), as n —> 0, 

f T (u K (s), C( S )e 2 — x ) ds -> / T (u( S ), r' m ( S )e 2 ™) ds, 
Jo Jo 

v C \(u K ( S ),T m ( S )e 2mm - x ))ds -> v / T ((u(.s),r m ( S )e 2 " m - x ))d S , 
Jo Jo 

[ T (6 K (s)e 2 , r m ( S )e 2 " m x ) ds -> f (9{s)e 2 , r m ( S )e 2Tlm x ) ds, 
Jo Jo 

V(«). e 2 ™)*^) ds H- f T (9(s), e 2 — ) x ' m (s) ds, 

o 

T 

27rim-x , 



((0 K ( S ),e~* Xm (s)))d S 



o 



< c^{^\\e K \\ L 2 ([0tTlm) ) < cko^^o 



It remains to show the convergence of the remaining non-linear terms. Let 

'("0 : =E/ («X,r m ( s )9 3e 2 " mx )d s -V (^u,r m (s)a,- 

3 = 1 J ° 3=1 J ° 

T pT 



e 2 " m ' x ) 



J(«):= / K( S )0 K ( S ), Xm ( S )Ve^ mx )d S - / (u(.s)<?(.s), Xm (s)Ve^ m ' x )ds 
Jo Jo 

The convergence I(k) —> as k —> is standard in the theory of the Navier-Stokcs 

equations, thanks to (3.5a) and (3.5b) (see, e.g., [16, 52]). To show J(k) — > as 

i% — ^ 0, we write J(k) = Ji(k) + J2(«0, the definitions of which are given below. Wc 
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have 



((u re (s) - u(s))9 K (s), Ve~ x ) Xm (s) ds -+ 



as k — > 0, since u K > u strongly in L 2 ([0, T],H) and 8 K — >■ 8 weakly in L 2 ([0, T],if). 
For J 2 , we have 

J a («) := f (u(s)(^(s) - 0(a)), Ve 2 ™ mx ) Xm («) <k -> 
Jo 

thanks to the weak convergence in (3.5d) and the fact that u e L 2 ([0, T},H). Thus, 
J(k) = Jx(k) + J 2 (k) —> 0. Hence, sending k — > 0, we see that u and 8 satisfy (3.1). 

Step 3: Show that u e C([0,T], #), t/wrf <9 G C«,([0, T], L 2 ), and that in fact 
0eC([O,T],i 2 ) . 

The uniform bound with respect to k on the time derivative of u K given in 
Theorem 3.2 (ii) allows us to pass to an additional subsequence if necessary to 
find that ^ G L 2 {[0,T},V). Since u g L 2 ([0,T},V) and ^ G L 2 ([0, T],V'), 
following the standard theory of NSE, (see, e.g. Theorem 7.2 of [44]) we obtain 
that u e C([0,T],H). 

Next, we would like to show that 8 € C w ([0, T],L 2 ). This can be proven without 
difficulty using standard arguments. For completeness and for use in the later 
section we present the proof here. We follow similar arguments as in [37, 52]. 
We start by showing that the sequence of solutions {8 K } (as k — > 0) is relatively 
compact in C w ([0,T],H). By the Arzela-Ascoli theorem, it suffices to show that 
(a) {8 K (t)} is a relatively compact set in the weak topology of L 2 ([0,T],T 2 ) for 
a.e t > and (b) for every £ L 2 (T 2 ) the sequence {(8 K , 4>)} is equicontinuous in 
C([0, T]). Condition (a) follows from the uniform boundedness of 8 K (t) in L 2 (T 2 ) 
for a.e. t > 0, as stated in Theorem 3.2 part (i). Next, we show that condition (b) 
is satisfied. Following classical arguments, we start by assuming that <f> is smooth, 
for example we can assume that <fi is a trigonometric polynomial. We have 

|(0 K (t 2 ),0)-(0 K (*i),0)| 



< 



(3.6) 



{{6 K {t),<£)) dt 
\8 K (t)\\\W(j>\ dt 



< K 



(V(u K K )(i),0) dt 

L 
j 

(u K (t)d K (t),Vcf>) dt 
1/2 



|V0||oo|*2 — *l| 



t 2 x V* 

4 



\u K (t)\\ldt 



K (t)\ 2 dt 



1/2 



Without loss of generality assume < k < 1 and use (2.14), one then obtains 

ft, \ V2 

\(8 K (t 2 ), <fi) — (8 K (ti), 4>)\ < C\t 2 — t^ 1 / 2 ( K I 

(3.7) 



\(8 K (t 2 ),<t>)-(6 K (ti)A)\ <C\t 2 -h\V 2 U I 

+ C\t 2 -t 1 \ 1 / 4 I |u K (i)| 2 ||u K (i)|| 2 di. 



ti 
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From Theorem 3.2 part (i), since u K is uniformly bounded with respect to k in 

L°°{[0 7 T],H) n L 2 ([0,T],V) and n [ \\9 K (t)\\ 2 dt < K , with K independent of 

Jo 

K, we have that the set {(9 K , 4>)} is equicontinuous in C([0, T]). We now extend this 
result for all test functions <j> in L 2 (Y 2 ) using a simple density argument of trigono- 
metric polynomials in L 2 (T 2 ). Let e > 0. We choose a trigonometric polynomial 
4> e such that \<p — e | < 3| < y_[_ 1 - Then, we have 
(3.8) 

\(6 K {t 2 ),<f>)-{0 K (ti),<t>)\ = |(^(i 2 )-^(ti),0-^) + (^(t 2 )-^(ii),0e)| 

< 10 - 4> 6 \ (\e K (t 2 )\ + 1^(^)1) + \(e K (t 2 ) - e K (h),cf> e )\. 

From the uniform L°° ([0,T], L 2 ) bound of 9 K , with respect to k, we conclude 
that the first term on the right-hand side of (3.8) is less than |e. Choosing | *2 — ^i| 
small enough in (3.7) we can make the second term on the right-hand side of (3.8) 
to be less than e/3. Thus, the whole expression can be made less than e. This 
completes the proof that 8 G C w ([0, T],L 2 ). Finally, as pointed out by the authors 
[17], since 9 is transported by the div-free velocity field u G L 2 ([0 7 T], V), we get in 
addition that 9 G C([0,T],L 2 ), (see, e.g. [20] ). 

From these results, standard arguments from the theory of the Navicr-Stokcs 
equations (see, e.g., [16, 44, 52]) now show that the initial conditions are satisfied 
in the sense of Definition 3.3. 

Step 4: Show that if 9 G L°° then 9 G L°°([0, T], L°°). 

Here we will use E. Hopf and G. Stampacchia technique which are very similar 
to those used in [24] (see also [32, 51]), but we give the details here for the sake 
of completeness. For any function / G H 1 , we use the standard notation /+ := 
max{/, 0}. It is a standard exercise to show that if / G H 1 , then /+ G H 1 . Let 
(u K ,9 K ) be a solution of (1.1), as given in Theorem 3.2. Let us denote K := 
9 K — ||#o||l°°- Notice that 9 K satisfies the evolution equation (1.1c) with 8 replaced 
by K and u replaced by u K . Thus we have (6 K ) + G £ 2 ([0, T], H 1 ). Taking the 
action of (1.1c) with (Q K ) + yields 

lj t W(®-) + Wh = -«/ 2 |v(e K )+| 2 dx + ^u K e K • ve+dx 

= -k[ |V(6 K )+| 2 dx+i / u K -V(6 + ) 2 dx 

= -k [ |V(6 K )+| 2 dx<0, 
Jt 2 

thanks to (1.1b). Thus \\(Q K )+(t)\\ L 2 < \\(Q K ) + (0)\\ L ' = which implies that 
(# K (x, t) — ||#o||l=°) + < a.e. Similarly, one can show that (H^oHi 00 —8 K (x, t)) + >0 
a.e. It now follows that H^IIl^^t],^) < H^olU 00 for all n > 0. Thus, we 
have 8 K is bounded uniformly with respect to k in L°°([0, T], L°°). Therefore, 
it follows from the Banach-Alaoglu Theorem, that there exists a subsequence of 
the previous subsequence which we also denote as 8 K converging in weak-* topol- 
ogy of L°°([0, T], L°°) to 8 and satisfies the following bounds: ||0||l°°([o.t].l°°) ^ 

lim i^ f ll^«IU oo ([0 1 T],£°°) < 11^0 |U« < 00. 

The equivalence of (3.1) to the functional form (3.3) follows from the standard 
argument of NSE (see, e.g., [52]). 

□ 
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Theorem 3.6 (Existence of strong solutions). Let uo <E V and 6q € L 2 . Then 
there exists a strong solution of P® . Furthermore, the functional equation (3.3a) 
now holds in L 2 ([0,T], H), and (3.3b) now holds in £ 2 ([0, T], H~ v ). 

Proof. Since Uo € V and do G L 2 , we have by Theorem 3.2 that there exists a solu- 
tion (u«,0«) of (1.1) with u K £ L°°([0,T],V) n i 2 ([0, T],T>(A)), and furthermore, 
€ L 2 ([0,T], H). In order to show that in fact the bounds on higher-order 
norms arc independent of k (as stated in Theorem 3.2 part (h)), let us take the 
inner product of (1.2a) with Au K . Using the Lions-Magenes lemma (see, e.g., [52]) 
to show that (^,Au) = |^j|u|| 2 and the fact that (B(u K , u K ),Au K ) = due to 
the periodic boundary conditions, we have 



1 d .I 
2^ K 



1 



-v\Au K \' = 

Subtracting -||Au K | 2 and using Gronwall's inequality yields 



K e 2 ,Au K ) < \9 \\Au K \ < —\9 \ 



^u K | 2 . 



(3.9) || 



u,. 



1 



1 



\Au K \ z ds < ||u |r + -\8 a \ z t < ||u f + := K 3 . 

Thus u K is bounded in L°°([0,T],V) L 2 ([0, T],T>(A)) independently of k. Fur- 
thermore, 



du K 



dt 



= sup (B(u K ,u K ),w) 

|w|=l 

< sup |Au K |||u K |||w| -I 

|w|=l 



- v sup (Au K ,w) 

|w|=l 

v sup |j4u k ||w| 

|w| = l 



- sup (tf K e 2 ,w) 

|w|=l 

sup |0 K ||w| 

|w|=l 



<if 3 |^u K | + ^u K | + |0 o 



Thus, is bounded in L 2 ([0,T], H) independently of k due to (3.9). We also 

have, 

dd K 



dt 



= sup |<£(u«,0 K ),w) 

if-i || w|| =1 

< sup ||u w ||i«.|0 re |||w 

II W II — 1 



« sup |(V0 K ,Vw) 

llwll — 1 



^ sup I 

llwll — 1 



w 



< ||u k || h2 |0 o | + v^I|0kII, 

where we have used here the assumption that < k < 1. Hence, from Theorem 
3.2, we have that ^jf is bounded in £ 2 ([0, T], H^ 1 ) independently of n. The above 
estimates allow us to use the Banach-Alaoglu Theorem and the Aubin Compactness 
Theorem (see, e.g., [16, 52]) , as K — > 0, to extract a further subsequence (extracted 
from the sequences in (3.5a), (3.5b) and (3.5d), and which we still label with a 
subscript K, such that 

(3.10a) 

(3.10b) 

(3.10c) 



(3.10d) 





■> u 


u« - 


A u 


du K 


du 


dt 


~dt 


de* 


a de 


dt 





weakly in L 2 ([0,T],H), 
weakly in L 2 ([0, T], i/" 1 ), 



where the limit u and are the same elements as in (3.5), by the uniqueness of 
limits since the current topology are stronger than those in (3.5). Furthermore, 
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since u G L 2 ([0,T],T>(A)) and ^ G L 2 ([0, T],H), following the standard theory of 
NSE, (sec, e.g. Theorem 7.2 of [44]) we obtain that u 6 C([0,T], V). Thus we have 
shown the existence of a strong solution as defined in Definition 3.4. □ 

In the next theorem we will show the uniqueness of strong solutions. We note 
that in the work of [18], global well-posedness in the case of whole plane R 2 was 
established with initial data Uo and Oq both only in L 2 . This optimal global well- 
posedness result was established using elegantly a priori estimates in Besov spaces 
for the heat equation and the transport equation. Here we give an alternate proof 
which uses more elementary techniques but we require stronger initial data for the 
velocity field. This will allow us to fix some basic ideas that we will use to get 
the optimal global well-posedness results for the anisotropic Boussinesq equations 
which we will present in the next section. 

Theorem 3.7 (Uniqueness of Strong Solutions of P® ). Let T > 0. Suppose 
0q G L 2 and Uq G V. Then there exists a unique strong solution (u,(9) to P® . 

Proof. The existence of solutions satisfying the hypothesis is already given by The- 
orem 3.6. It remains to show the uniqueness. Let (ue,6e) be two strong solutions, 
£=1,2, and define & := A" 1 ^ such that L 2 & dx = on the interval [0, T]. Write 

u := Ui — u 2 , and £ := £i — £ 2 . These quantities satisfy the functional equations 

(3.11a) — + vAu + B(ui,u) + B(u,u 2 ) = PjA^e 2 ) in L 2 ([0,T],H) and 
dt 

(3.11b) ^+£(u,A£ 1 )+£(u 2 ,A£) = in L 2 ([0, T], H' 1 ). 
dt 

Taking the inner product in H of (3.11a) with u , and taking the action in H^ 1 of 
(3.11b) on £ G L 2 ([0,T],# 2 ), wc obtain, thanks to Lemmas 2.1 and 2.2, 

(3.12a) ~\u\ 2 + u\\u\\ 2 = -(B(u, Ul ), 5) + (A£e 2 ,u), 

(3.12b) ~U\\ 2 = (uA£ lf V£) - (u 2 A£,V£). 

In (3.12b) we used the Lions-Magenes Lemma (see, e.g., [52]) to obtain ^^||£|| 2 = 

(^i'V' L0t K = maX£=1 < 2 {II u cIIl~([o,t],v), ||^||l°°([o,t],l2)}- From equation 
(3.12a), (2.8) and since Ui G L°°([0,T], V) we have 

illu^ + HISf^ClulllullKII + liailull 

(3-13) <£| 5 |» + £||u||' + |;||£|| 2 + ^||Sf. 

Next, let e > be given such that e <C 1. For the equation (3.12b), we integrate 
by parts and use (2.13) to find 
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Id 

2 dt [ 



--ii^ii 2 = (sAa,ve)+^(^u 2 ,v^ 

< ||u|| L »|A6|||eil + l|Vu a || iV .||f||||VfH 



L2/(l-») 



1— e || tile 

H 2 



^^(HuHe-^ + IAule- 1 /^ 4 ) ||£|| + C||Vu 2 || L2/c ||£|| ||£|M|£| 
(3.14) < ^||5|| 2 + C (^^ 1/2 + l) lieil 2 + \Au\ 2 e' 2 ^ 1/4 

+ ^ e - 1 /2|^ U2 |||e|| 2 - e , 

where we have used (2.2), (2.17), and the interpolation inequality || V£|| L 2/(i-«) < 
C||£ll 1_e ||£llij2, noting that C is independent of e. 

Next, we will use the fact that, £(0) = and u(0) = 0, and also that ||£(i)|| and 
|u(t)| are continuous in time and thus there exist a r > such that ||£(i)|| < 1 
and |u(f)| < 1 for all t e [0,t]. Let t* = sup{r € (0,T] : |u(i)| < 1 and ||£(t)|| < 
1 for all t G [0,t)}. Adding (3.13) and (3.14) and rearranging, we have on [0,f*], 

1 d 

lull 



(|uj 2 + ||£| 



2dtV J 2 

<kJi+ (|5| 2 + ||el 2 ) + \Aufe-^ 1/4 + K<e-V*\Au 2 \M\\< 



(3-15) <kJi + -^ + §^\Au 2 \\ (|S| 2 + Hell 2 ) 1 " + |Au| 



2 p -2/e 



1/4 



Let r\ > 0, be arbitrary and let z := |u| 2 + ||£|| 2 + r\. Dividing (3.15) by z 1 e , we 
find 



< ff„ (l + -L + ^\Au 2 \\ + (vy-^Aufe- 

since z>rj. Integrating over [0,t], for i £ (0, t*], we find 

/ /p \ ij 

(3.16) z(t) < 



1/4 



/^ /e f eT + e 1 / 2 ^ + K'e 1 ' 2 £ \Au 2 {s)\ ds^ 



+ 6 i/«( 7? )i-i/e e -a«-'-| I |^( S )| 2 d S +r/ 



Sending r? — > 0, we obtain 



(3.17) |S(i)| 2 + ||£(t)|| 2 < K x J e (^eT + e x ' 2 T + cK'e 1 ' 2 £ \Au 2 {s)\ 2 ds^ 

for t e [0,i*]. Taking the limit of (3.17), as e -> 0, we find that U(t)\\ = and 
|u(i)| = on [0,t*]. In particular, |u(t*)| 2 = U{t*)\\ 2 = < 1. Therefore, from the 
continuity of |u(i)| 2 and ||£(^)l| 2 and the definition of t* , we conclude that t* = T, 
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otherwise we a contradiction to the definition of t* . Hence, u(t) = and =0 
for allte [0,T]. □ 

4. Global Well-posedness Results for the non-diffusive Boussinesq 
Equations with Horizontal Viscosity (P° ) 

We now consider the Boussinesq equations with anisotropic viscosity as given in 
(1.3) (P° x ). We will establish here global well-posedness results under some not 
too restricted initial conditions. In the first part of this section we will first define 
what we mean by weak solution to system (1.3) and then show its existence. Then, 
under some additional requirements on initial data, we can show uniqueness. To 
set additional notation, we denote the vorticity u> := d\u 2 — fou 1 , which satisfies 
the following equation 

(4.1) d t LJ + V • (wu) - vd\u = di6. 

The best global well-posedness result we are aware of for problem (1.3) in the 
case of the whole plane R 2 is stated in following theorem, established in [17]. 

Theorem 4.1 (Danchin and Paicu,[17]). Let il = R 2 . Suppose 9q G L 2 L°° , and 
Uq G V with ujq G \~L. Then system (1.3) admits a global solution (u, ff) such that 
9 G C B ([0,oo);L 2 ) n C w ([0,oo);L°°) n L°°([0, oo), L°°) and u G C w ([0, oo); H 1 ), 
u e 2 G ^([O.oo);^ 2 ), u G L£ c ([0, oo), VI), Vu G L 2 OC ([0, oo), VI). If in 
addition 9 a G H s for some s G (0,1], then 9 G C([0, oo); H s ~ e ) for all e > 0. 
Finally, if s > 1/2, then the solution is unique. 

In the present work, we improve the above result by weakening the requirements 
on the initial data needed for the uniqueness portion of the theorem. To begin 
with, we weaken the notion of solution by making the following definition. 

Definition 4.2 (Weak Solutions for the Anisotropic Case). Let T > 0. Let 9 a G L 2 , 
ujq = V- 1 • u G L' 2 . We say that (u, 0) is a weak solution to (1.3) on the interval 
[0,T] if u e i oo ([0,T];L 2 )nC„,([0,T];i 2 ) and 9 G L°°([0, T}; L 2 ) n C w ([0, T}; L 2 ), 
u 2 G L 2 ([0,T],iJ 2 ), ^ G LHP,? 1 ],!/'), f € L\%T\,H- 2 ) and also (u,6») satisfies 
(1.3) in the weak sense; that is, for any $, <p, chosen as in (3.2), it holds that 

T i-T 2 t 

(u(a), *'(«)) ds + i/ / (5m(s),fli*(a))ds + V / («*u,fl!,-$) ds 
Jo J=1 Jo 

(4.2a) =(uo,$(0))+ / (0(.s)e 2 ,$(s))ds 
Jo 




(4.2b) -/ (6(s),<p'(s))ds+ [ (6u,Vip)ds = (0 ,v(0)), 
Jo Jo 

where ' = 4-- 

as 

Remark 4.3. Again following standard arguments as in the theory of NSE [52] one 
can show that the above system is equivalent to the functional form 

(4.3a) ^ + vd?u + B(u,u) = PJ9e 2 ) in L 2 ([0, T],V') and 

dt 

d9 

(4.3b) — +B(u,6) = in L 2 ([0, T], iJ" 2 ). 
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We now state and prove our main results for the system (1.3) (P® ). The 
global existence and regularity results will be stated in the theorem below and the 
uniqueness theorem will follow. 

Theorem 4.4 (Global Existence and Regularity). Let T > be given. Let Oq G L 2 
and luq G L 2 . Then, the following hold: 

(1) There exists a weak solution to (1.3) (Pj) ) in the sense of Definition 4-2. 

(2) If U)q G LP, and 9q G L p , with p G [2, oo) fixed, then this weak solution 
satisfies uj G L°°([0,T],L p ) and 9 G L°°([0,T],L p ). 

(3) Furthermore, if ujq G \fh and 9q G L°° , then there exists a solution to G 
L°° ([0, T] , sfl)r\C w {% T],L 2 ),|eI 2 ([0,T], V) and 9 G L°° ([0, T] , L°°)n 
C([0,T],w*-L°°) (where w*-L°° denotes the weak-* topology on L°° ) with 
f t Gl oo ([0,T],F- 1 ). 

Proof. The outline of our proof is as follows. We begin by generating approxi- 
mate sequence of solutions (u^ n \ 9^ n ') to P® by adding artificial vertical viscosity 

v^f 1 > 0, artificial diffusion > 0, where K^ n \ Vy — > as n — > oo, and also by 
smoothing the initial data. Global existence of solutions to the fully viscous system 
P" K , given smoothed initial condition is guaranteed (see, Theorem 3.2 part (hi)). 
Next, we establish uniform bounds, for the relevant norms of the approximate se- 
quence of solutions which are independent of n using basic energy estimates. We 
then employ the Aubin Compactness Theorem (see, e.g., [16, 52]) to show that 
the sequence of approximate solutions has a subsequence converging in appropriate 
function spaces. This limit will serve as a candidate weak solution. We then show 
that one can pass to the limit to show that the candidate functions satisfy the weak 
formulation (4.2). Then we establish some regularity results. 

Step 1: Generating solutions to the regularized system given smoothed initial 
data. 

Let v x > be fixed and let be a sequence of positive numbers, con- 

verging to zero. In fact, we can also assume that both < v x and < v x . 
Let (Uq™ , #o ) is a sequence of smooth initial data such that Ug™^ — > Uo in V and 
9^ —t 9q in L 2 , chosen in such a way that for each n G N, ||u ™' ) | < |uo|| + ^ 
and \9q \ < \6q\ + — . Notice, since Ug is smooth it follows that V -1 ■ Ug = oj^ 

and so Wg are smooth functions bounded in L 2 . From Theorem 3.2 part (iii), 
by slightly modifying the proof of this result to account for values of the viscosity 
which differ in the horizontal and vertical directions, we have that for each n, there 
exist (u^ n \9^) satisfying the following equations: 
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(4.4a) 

- f {v^\s),^'{s))ds + v x I (d 1 uW(s),d 1 $(s))ds 
Jo Jo 

+4 n) / (d 2 u^(s),d 2 ^(s))ds + V f (u^ n) u^\d^)ds 
Jo j =l Jo 

= ( U ("\$(0))+ [ T (e^( s )e 2 ^( s ))ds 
Jo 

(4.4b) 

-/ (6^(s),tp'(s))ds + n<< n) [ (W n \s),V$(s))ds+ [ (8^u( n \Vtp)ds 
Jo Jo Jo 

(4.4c) =(<#°,¥>(0)). 

Step 2: A priori estimates and using compactness arguments to prove convergence 
of a subsequence. 

We next establish a priori estimates on (vS n ',6^ n ') uniformly in n (independent 
of v y and From the above smoothness properties of (u^™', £>(")), we can now 

derive a priori estimates using basic energy estimates in which the derivatives and 
integrations are well defined. First, one can obtain, because div u( n '=0, that 

(4.5) \0 {n) (t)\<\0l ) n] \<\e o \ + -, 

n 

and 

|u< n >(t)| 3 +2i/ x f |a lU (")(r)| 2 dr + 2^") f \d 2 u^ (t)\ 2 dr 
Jo Jo 

<(\u \ + - + t(\e Q \ + -)) 2 . 

n n 

The calculations above are justified by replacing the test functions by 9^ and uy 1 ' 
in (4.4) and then integrating by parts. 

Using the evolution equation of the vorticity, namely the equation 

(4.6) d t J n) + u (n) • Vw (n) - p x dlJ n) - i/W^«W = 6 x n) , 
we also have 

<^ia lW («)|2 + _L|^)|2_ 

Integrating this gives 

(4.7) \u> (n) \ 2 + v x [ \dnj {n) \ 2 dT + 2is y n *> f \d 2 ^ (n) \ 2 dT 

Jo Jo 

1\ 2 t / , l^ 2 



(4.8) <(M + ") +2^W + n 
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which implies that is uniformly bounded in L°°([0, T], L 2 ) with respect to 
n, and therefore uW is uniformly bounded in L°°([0,T],V) with respect to n. 
Furthermore, (4.7) shows that d\ ojW is uniformly bounded in L 2 ([0,T],L 2 ) with 
respect to n. We also observe that 

Therefore, Am 2, W is uniformly bounded in L 2 ([0, T], L 2 ), so that u 2 '( n > is uni- 
formly bounded in L 2 ([0, T],H 2 ) by elliptic regularity, and thus Vu 2 >(") is uniformly 
bounded in L 2 ([0,T], H 1 ), all with respect to n. Next we derive uniform bounds 
on the derivatives ( d "' t - ) rag M- Note that 

J, 



Thus, 



<7f 



< sup 

||w|| f v 2 =l 



i/W sup 



(4.9) 



sup 

l|w|U 2 =l 



,2=1 



i/W sup 



sup 

I wll v,2 =1 



sup 

I|W|| A 2=1 

sup 

l|w|| A 2=l 

sup 

||w|U 2 =l 



ai6>( n) ,w 

6»^,5iw 



< |w(")||u(")| 1/2 ||u (ll) H 1/2 + i>x|w< n) | + i/ x |wW| + |6>W|, 

Since each of the terms on the right-hand side of the inequality above is bounded 
independently of n, we deduce by the Calderon-Zygmund elliptic estimate (2.3) 
that <9 t u(") is bounded in L°°([0,T],V') independently of n. Similarly, one can 
show easily that 

(4.10) < |^")||uW| 1 / 2 ||uW|| 1 /2 ) 

H- 2 



dt 



which implies also that at> dt ' is bounded in L°°([0, T],H 2 ) independently of n. To 
summarize, we have from the above results that 



(4.11a) 




is 


bounded 


in 


L°°([0 t T\,L 2 ), 


(4.11b) 


(uW)„ eN 


is 


bounded 


in 


L°°Q0,T\,V), 


(4.11c) 




is 


bounded 


in 


L 2 ([0,T],# 2 ) 


(4.11d) 


V dt Jnen 


is 


bounded 


in 


L°°({0,T],V), 


(4.11c) 


ue^\ 

V dt Jnen 


is 


bounded 


in 


L°°([0,T},H- 2 



Using Banach-Alaoglu and Aubin Compactness theorems (see, e.g., [16, 52]), the 
uniform bounds with respect to n as stated in (4.11) implies that one can extract 
a further subsequence (which we relabel with the index n if necessary) such that 
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(4.12a) 6 {n) weakly in L 2 ([0,T],L 2 ) and weak-* in L°°([0, T], L 2 ). 

(4.12b) u (n) -> u strongly in L 2 ([0, T},H), 

(4.12c) u (n) u weakly in L 2 ([0, T], F) and weak-* in L°°([0, T], V), 

(4.12d) w 2 ' (Tl) u 2 ^ (,l) weakly inL 2 ([0,T], H 2 ), 

(4.12e) — — weakly in L ([0, T\,V') and weak-* in L°°([0,T],V'), 

\Jb\j \Jb\j 

(4.12f) — >■ — weakly in L 2 ([0, T], #" 2 ) and weak-* in £°°([0, T},H~ 2 ). 

Step 3: Pass to the limit in the system. 

It remains to show that (4.12) is enough to pass to the limit in (4.4) to show that 
(u, 9) satisfies (4.2). To do this, in accordance with Remark 3.1 and Definition 4.2, 
we only consider test functions of the form (3.2), which we note is sufficient for 
showing that (u, 0) satisfies (4.2). For the linear terms in (4.4), we have, by the 
weak convergence in (4.12c) and (4.12a), asn->oo (that is, K^ n \vy — > 0), 

[ T (uM( S ),T' m ( S )e 2 —)ds -+ / T (u(.s),C( S ) e 2 — x )ds, 
Jo Jo 

Vx / T (a lU (")( s ),r m ( s )9 1 e 27rim - x )d s -> v x / T (o lU ( s ),r m ( s )a ie 2 " m - x )d s , 

Jo Jo 

[ T (8^(s)e 2 ,T m (s)e 2 ™ m -*)ds^ f ^(.s)e 2 , r m ( S ) e 2 ™ m x ) da, 
Jo Jo 

(0^( S ),e 2 —-) X ' m (s)ds -+ [ T (0( S ),e 2 ^) X ' m ( S )d S , 

Jo 

K (n) f 1 (d2U^( S ),r m ( S )d 2 e 2mm -*)ds -> 0, 
Jo 

((6 K (s),e 2 ™ m -* Xm (s)))ds 



< CKqVkW -> 0. 

It remains to show the convergence of the remaining non-linear terms. Let 
2 „t 2 T 

J(n) := V / (« J ''WuW > r m (s)fl!je 2 ' r4m - x ) ds-J^ (tfu.r^s)^ 2 ™' 1 ) 

J(n):= / (u(")( s )0(")( s ),x m ( S )Ve 2 " m ' x )d S - / (u(*)0(*), x m (s)Ve 2 ™ m ' x )<fe 



22 



ADAM LARIOS, EVELYN LUNASIN, AND EDRISS S. TITI 



To show I(n) — > as n — > oo, we write I(n) = Ii(n) + I^in), the definitions of 
which are given below. We have 



\h(n)\ :- 



< 



f T ((u j ' (n) (s) - u-'»)u(")( S ),^e 2 " m - x ) Xm ( S )d S 

7=1 J ° 



|u(")( S )-u( S )||u(' l )( S )||Ve 2 " m ' x x m ( S )|d.s 



|| U W - U|| r .2 „_ HU^II L ¥ H* || Ve 2 " m - X Xm|U|L r 0, 



asn-> oo, since u^ 1 ' — > u strongly in L 2 ([Q, T], H) and 

u («) i s 

uniformly bounded 

in L°°([0,T],V) and hence in L°°([0, T],H) . Similarly, for 7 2 , we have that as 

n->oo 

2 -T 

/ 2 (n) := £ / (^'(s)(u (n) ( S ) - u( S )),a, e 2m;m - x ) Xm (s)ds -> 0. 

3=1 ^° 

To show J(n) — > as n — ?> oo, we write J(n) = J\{n) + J2(n). We have 
h{n) := f T ((uW( S )-u( S ))^»(s) I Ve 2, " m ' x )Xm( S )rf^0, 



oo, since u 



u strongly in L 2 ([0, T], iJ) and 9^ -> 9 weakly in 



L 2 ([0,T],H). For J 2 , we have 



J 2 (n) := / (u(s)(0W(«) - 0(a)), Ve 2 ™) Xm (s) 



ds -> 0, 



by the weak convergence in (4.12a) and the fact that u E L 2 ([0, T],H). This estab- 
lishes the existence of weak solution to the system P° when uo G H 1 and 9q E L 2 . 

Step 4: STiow that oj E C w ([0, T];L 2 ). 

By the Arzela-Ascoli theorem, it suffices to show that (a) {ui^} is a relatively 
weakly compact set in L 2 (T 2 ) for a.e t > and (b) for every <f> E L 2 (T 2 ) the 
sequence {(u/™),0)} is cquicontinuous in C([0,T]). Condition (a) follows from the 
uniform boundedness of wW in L 2 (T 2 ) for a.e. £ > given in (4.7). Next, we show 
that condition (b) is satisfied. We follow similar argument as in Step 3 of Section 
3 equation (3.6), where, we start by assuming that is a trigonometric polynomial 
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to obtain, 

|( W W(ta)^)-(a;W(*i)^)l 

< \ Vx [ V« (n) (*)>3i0 dt\ + \u v [ ' \d 2 ^ n \t),d 2 cp) dt\ 

Jti Jti 

+ 1 f(uW.v^, w W) dt\ + 1 f\e^\d x <t>) dt\ 
Jti Jti 

< v x t \d^\\d x 4>\ dt + v x [* \u^\\dU\ dt 

Jtl Jtl 

+ [' |u(")|| w (™)| di + / 2 \0^\\V(/>\ dt 

Jti Jti 

< |WU|ta - ti| 1/2 ^ f 2 |9l^^| 2 + |fi^U|t a - tilllw^Hi-ij 

+ ||V0|U|t a - h\ (||uW|| i - ?£ ;||wW|| i? . £2 + ll^lU-ij) , 

where recall we have assumed without loss of generality that v y < v x . From 
the uniform boundedness of u)( n > (4.7) and 9^ (4.5), the right-hand side can be 
made small when \t 2 — t\\ is small enough. Thus we have that the set {(a/"', <fr)} is 
equicontinuous in C([0, T]). Then one can extend this result for all test functions 
(f) in L 2 (T 2 ) using a simple density argument as before. This completes the proof 
of part (1) of Theorem 4.4. 

Step 5 Proof of part (2) of Theorem 4-4- 

We choose a sequence of smooth initial data <Jq — > loq and similarly 9^ — > #o 
in every LP with p > 2 chosen in such a way that for each n <G N, ||wq || p < 
||wollp + ^ an d ||^o lip — ll^ollp + h- From Theorem 3.2, we obtain for each n, a 
solution %v- n > G H 3 which then gives us a/™-* g H 2 which is a topological algebra, 
hence |u>) |?>- 2 u>) g ff 2 . We take the inner product of (4.6) with \uj^ \v~ 2 u^. 
Integrating by parts, we have 

-4ll"> (n) IIS + *'*(P-l) / \dxu {n) \ 2 \u>^\v- 2 dx + vt n \p-l) [ \d 2 ^\ 2 \u\v- 2 dx 
pat p J T 2 y J T 2 

<(p-l) / \9 < -^\\d 1 J n) \\uj {n) \ p ~ 2 dx 
Jt 2 

<v x (p-l) I \d 1 J n) \ 2 \^ n) \ p - 2 dx + ^—^- f \9 { ^\ 2 \J^\P- 2 dx 
Jt 2 Jj2 

<u x (p-i) [ |5 lW (")| 2 i^")r 2 dx+^ii^)|| 2 i| W (")||p- 2 . 

Jt 2 4v x 

Therefore, we have 

^M" , ii?^ii» w i&i" w ir 2 ^(ii*ii»+;) 2 ii" w r 2 - 

That is, 

2 
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Integrating in time, we have 

(4.13) \\^ n \t)\\l < \U n) \\l + ^ (\\0o\\r + 



< 



flkollp H 




_v-U 








.i)'. 




INI P H 


-i) 


n J 


2u x \ 


n / 







_p-U 








.i)'. 




IN| P H 


-i) 


n J 


2v x \ 


n / 



That is, w' n ) is uniformly bounded in L°°([0, T], L p ) for each p G [2, oo), inde- 
pendent of n. It follows from the Banach-Alaoglu Theorem and diagonalization 
process, that there exists a further subsequence which we also denote as a/"' con- 
verging weak-* in L°°([0, T], LP) to some limit which we denote as w and this limit 
also enjoys the limit of the upper bound, that is 

/ \ 1 / \ 2 

(4-14) \\u\\l< 
This implies that w G i°°([0, for all p G [2, oo). Similarly we find that 

(4-15) ||^ ) (*)ll P <ll^ n) || P <ll^o|| P + ^ 

which implies that 0(") converges weak-* in L°°([0, T];L P ) to 8 £ L°°([0, T];L P ) 
for allpG [2,oo), and ||0||i,.*.([o,t] 1 l:p) < ll#o|| P - 

Step 6 Proof of part (3) of Theorem 4-4- 

To prove Theorem 4.4 part (3) we divide both sides of (4.14) by p— 1 and then 
taking the supremum over all p > 2 of both sides, we get that w G L°°([0, T],y/~L) 
provided that wo & \/L and 8q G Next, we want to show that G C([0, T]; w*-L° 
We will use the Arzela-Ascoli theorem as in Step 4. Notice that if 8q G then 

(4.15) holds uniformly for all p G [2, oo) and hence 

(4.16) l|0 (n) (*)IU < ||0o||°o + -- 

n 

This implies that the sequence 8^ n \t) is a relatively compact set in the weak—* 
topology of L°°([0,T] x T 2 ). It suffices to show that the sequence {(6»(«), </>)} is 
equicontinuous in C([0,T]) for every (f> E L 1 . It follows automatically from the 
previous result and the density of L 2 (T 2 ) in L 1 (T 2 ) that 9 G C„,([0, T], L 2 ). Finally, 
we would like to show that § G L°°([0, T], H" 1 ) and hence § G L 2 ([0, T], ff- 1 ). 
Since w G L°°([0, T], vT), we have in particular that w G L°°{[0, T], L 3 ), and 
hence u G I^QO, T], W 1 ' 3 ) C L°°([0, T], L°°) by (2.3), (2.16), and the Sobolev 
Embedding Theorem. From equation (4.3b), using (2.11) and the fact that 8 G 
L°°([0,T],L 2 ), we obtain, 

(4.17) ^ = sup |(B(u,ff),«;>|<||u|| oo |fl|<ooa.etG[0 i r]. 

This completes the proof of part (3) of Theorem 4.4. □ 

Theorem 4.5 (Uniqueness for the Anisotropic Case). Let 8q G L°° , ljq G \f~L. 

Then, for every T > 0, there exists a unique solution w G L°°([0, T], \[~L) n 
C w ([0,T];L 2 ) and 8 G L°°([0, T], i°°) n C*([0, T]), w*-i°°) to (1.3). 
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Proof. Let T > arbitrarily large. The existence of solution on the interval [0, T] is 
established above, therefore it suffices to show uniqueness. We note that some very 
important a priori estimates that we need in the beginning of this proof were first 
elegantly derived in [17]. We recall those estimates that we have borrowed from 
[17]. We have derived them rigorously in the previous theorem and we derive them 
here again formally to make the proof of uniqueness self-contained. First, one may 
easily show that for any pg [2, oo], we have 

(4-18) \\6{t)\\ v < \\6 Q \\ p , 

so 6 G L°°([0, T],L P ), p G [2, oo]. Given that cj G V% and hence loq G L 2 , we have 

I| M 2 + v \% u \ = -{o&u) < v -\dM 2 + ^\o\ 2 . 

Integrating this gives 

\uj\ 2 + v f \8M 2 dr< | Wo | 2 + -|0 o | 2 . 
Jo v 

This implies that oj G L°°([0, T], L 2 ), and therefore u G L°°([0,T},V). Further- 
more, diui G L 2 ([0, T], L 2 ). Using the divergence free condition (1.3b), we observe 
that 

diu = d\u 2 - dxd 2 u x = d\u 2 + d 2 2 u 2 = Au 2 . 

Therefore, Au 2 G L 2 ([0,T},L 2 ), so that u 2 G L 2 ([0, T],H 2 ) by elliptic regularity, 
and thus Vu 2 G L 2 ([0, T\,H l ). By inequality (2.2), we have 

(4.19) ||W!| P < Cv^ - "T||Vu 2 ||^i. 

so that Vu 2 G L 2 ([0,T},VI). 

Next, we recall that we have global in time control over the ||a;||^. Taking the 
inner product of (4.1) with |a>| p-2 u; for some p > 2 and integrating by parts, and 
integrating in time, we have 

(4-20) \\ u (t)\\l <\\"o\\ 2 p + ^\\eo\\ 2 p t. 

This shows that u> G L°°([0, T], V~L). Using this, and the facts that diu 1 = 
—d2U 2 (by (1.3b)) and c^m 1 = d\u 2 — w, we have thanks to (4.19) that Vu 1 G 
L 2 ([0,T], y/Z). Combining this with (4.19) shows that 

(4.21) VuG L 2 ([0,T],\/L). 

We recall again that all the estimates above were first derived in [17] for the case 
where Q = M 2 . 

We are now ready to show that if (ui, 9\) and (112, 02) are two solutions to (4.2) 
on the interval [0, T] , with the same initial data (uo, #0) then they must be the equal. 
Define u := Ui u 2 , 9 := 61 — 9 2 , and & := A _1 #£, i = 1,2, and £ := £1 — £2. 
Based on Remark 4.3, these quantities satisfy the following functional equations. 

(4.22a) -£ + i/fr 1 u + fl(u,Ui) + fl(ua J u)=P ff (A£e 2 ) in £ 2 ([0,T],U') and 
at 

(4.22b) ^+£(u,A£ 1 ) + £!(u 2 ,A£) = in L 2 ([0, T], iT 1 ). 
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Taking the action of (4.22a) on u in L 2 {[0,T},V) and the action of (4.22b) in 
L 2 ([0,T],ff _1 ) with £ e L 2 ([0,T],H 2 ), thanks to the properties of the operator B 
in Lemma 2.1 and the operator B in Lemma 2.2 we obtain the following: 

2 

|5(i)| 2 + H|di5|| 2 = 5^(t?'ui,a,-u) + (A£e 2 ,u) 

i=i 

\j t \mf = "(SAfi, VO - (u 2 A£, v£), 
where again we have used Lions-Magenes Lemma (see, e.g., [52]) to get that u) = 
i^|u(i)| 2 and = hftU^W 2 - B y Lemma 2.1, we obtain 



ld_ 

2d? 



1 d 
2dr 



5| 2 + Hdi5| 2 < / |V Ul ||S| 2 dx+ (A£e 2 ,G) 



T 2 



and 



--U\\ 2 < 



T 2 



u • V£A£i dx 



u 2 • V£A£ dx 



Next, observe that, due to the divergence free condition, ei • diu = — e 2 • <9 2 u, we 
have 

|(Afe 2l u)| < / (\d^e 2 ■ &u| + |d 2 £e 2 • 9 2 S|) dx 



T 2 



|<9i£e 2 • 9m| + |92^ei • 9iu| dx 



< -IdxCI 2 + V -\e 2 ■ d,n\ 2 + ~\d 2 e + jlei • 9i5| 2 . 
z/4 i/4 



Combining the above estimates, we find 
Id 
2dV 



u\ z + v\d!u\ z < / ivuii isr dx + -neir + d^i 



< u 



|2/P 



T- 



ivuiiiar^x + ^i^ + ^isi 2 

^ 2 



^llVuxllpHSII^ISr^ + ^llef + ^xGF 

where we have used Holder's inequality. Similarly, by Lemma 2.2 
1 d 



2 dt 



u • V£A£i dx 



T 2 



|Vu 2 ||V£| 2 dx 



|2/p|r7<r|2-2/p 



< |u||V^|||Aa||oc + ||Vu 2 || p ||Var|V^ 

From the estimates above we can now adapt the well-known Yudovich argument for 
the 2D incompressible Euler equations (see, e.g., [27]) to complete the uniqueness 
proof. Let X 2 := |u(t)| 2 + ||£(£)|| 2 + ?/ 2 for some arbitrary r) > 0. Adding the above 
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2-2/p 



two inequalities and using Young's inequality gives, 

2dt 2 1 1 

<A^(|S| 2 + ||e|| 2 +r? 2 ) 

+ (||Vu 2 || p + HVuxllp) (||u||^ + HVeil^) (|S| 2 " 2/P + IVei 

< K V X 2 + C (||Vua|| p + HVuill,) (||u||^ + ||v£j|^) X 2 -*'*. 

Neglecting the term ||<9iu| 2 , dividing by X, and making the change of variables 
Y(t) = e~ Kvt X(t), we have after a simple calculation, 

Y < Ce- 2 *""* (||Vu 2 ||p + j|V Ul || p ) (||u||^ + ||VC||^ P ) F 1 " 2 ^. 

Integrating this equation and using the fact that e~ 2K " t ^ p < 1, we get that 

r r*1 / x lP/a 

K(*)< ^ 2/p + C/ -(llVuaWllp + IIVuiWyriluWII^ + llV^II^) d S 
Jo P v y . 

Letting 77 — > we discover that for all t G [0, T], 

|SW| 2 + ll^)ll 2 <fl|5|k^ + ||Ve||^L 



(4.23) • (c J X - (||Vu 3 («)|| p + ||Vm(s)|| p ) ehj . 

Thanks to the fact that A£ = 6 G L°°([0, T], L°°) C L°°([0, T], Z 4 ), we have by 
elliptic regularity that £ € L°°([0, T], W 2 ' 4 ), and thcrcforcjv 7 ^ G £°°([0, T], W 1 ' 4 ). 
Thus, by the Sobolev Embedding Theorem, we have V£ G L°°([Q, T], W 1 ' 4 ) C 
£°°([0,T],C ' 7 ), for some 7 G (0,1). Furthermore, 5 G L°°([0, T], \fZ) implies, 
for instance that u G L°°([0, T], W 1 ' 4 ) by the Calderon-Zygmund elliptic estimate 
(2.3). Using the Sobolev Embedding Theorem again, we have u G L°°([0, T], C° n ), 
for some 7 G (0,1). Therefore, the first factor on the right-hand side of (4.23) 
is bounded. Now, since Vu f G L 2 ([0, T], \fZ), £ = 1,2 by (4.21), we have by 
Cauchy-Schwarz 

" V "^" p ds < (t f T s»J VMS ^ ds) 1 ''. 

P \ JO P>2 P - 1 J 

f T ||VU £ (S)|| 2 

Let M e = / sup —-^ ds, £=1,2 and M = max{M 1; M 2 }. Thus, from the 

Jo P>2 P - 1 
above, for every fixed r G (0, T] we have 

(4.24) \u{t)\ 2 + ||£(t)|| 2 < K(2CMt) p / 2 , for all t G [0,r], 

where the constant C is the same constant which appears in (4.23) and K = 
(||u||l° ? i,~ + ||Vf||£«L»). Now choose r = r = min{T, ^j}, and consider 
(4.24) on [0,t ]. Taking the limit as p -> 00, we get that |u(t) | 2 + |j£(^)|| 2 < for 
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all t £ [0, tq]. Restarting the time at t = tq and noting the fact that 



t+T °\WMs)\\ Pds< L [ T snp \WM S )\\l ds 

P \ Jo P>2 P-l 



1/2 



we obtain from the analogue of (4.23) on [t ,T] that \u(t)\ 2 +\\£,{t)\\ 2 < K(2CMt ) p/2 
for all t £ [to, 2to]. Since we defined ro < 4 J A/ , we take the limit p — > oo and find 
that on the interval [ro,2ro], we also have that |u(t)| 2 + ||£(i)|| 2 < 0. We can 
continue this argument on the intervals [2to,3to], [3to,4to], . . . , and so on. Thus, 
we have |u(i)| 2 + ||^)|| 2 < for all t £ [0,T]. This implies that, \u(t)\ = and 
||f(t)|| = for alU £ [0,T\. □ 

5. Global Well-posedness Results for the Voigt-regularized Inviscid 
and Non-diffusive Boussinesq Equations (Poo) 

In this section, we investigate the problem P^q, a > 0, given by (1.5) (with v = 
K = 0) in 2D. We first establish global well-posedness results, and then investigate 
the behavior of solutions as a — > 0. In particular, we compare the limiting behavior 
to sufficiently regular solutions of the P§ problem. This leads to a new criterion 
for the blow-up of solutions to the Po problem. A similar criterion was given 
for the blow-up of the Surface Quasi-Geostrophic equations in [31], for the Eulcr 
equations in [34], and for the inviscid, resistive MHD equations in [33]. 

Definition 5.1. Let T > 0. Suppose Uo £ V and 9o £ L 2 . We say that (u,9) is a 
weak solution to the problem Pq on the interval [0, T] if for all test functions $, 
ip chosen as in (3.2), (u, 6) satisfies 

rT P T 2 

(u(s),&(s))ds - a 2 

.In 

3 = 

(5.1a) = (u ,$(0))+a 2 ((u ,$(0)))+ / (0(s)e 2 , *(«)) da, 



[ ((u(s) I * / («)))ds + V f (« J 'u, fl,-*) da 
Jo =1 Jo 



j (6(s),(p'(s))ds+ [ (6u,V(p)ds = (6 ,(p(0)). 



and furthermore, u £ C([0,T],V), !§ £ L°°([0, T],V), 9 £ L°°([0, T], L 2 ), 9 £ 
C w ([0,T],L 2 ) and § £ L°°([0, T],H~ 2 ). 

Remark 5.2. Following similar arguments as those for the NSE presented in [52] 
one can show that this definition is equivalent to the functional equation 

(5.2a) (I + a 2 A)— + P(u,u) = P a (9e 2 ) in L 2 ([0, T], V') and 

dt 

d9 

(5.2b) — + B(u,9)=0 in i 2 ([0, T],H~ 2 ). 

dt 

Theorem 5.3. Let Uo £ V, 9o £ L 2 . Then there exists a solution to Pqo' ^ n the 
sense of Definition 5.1. Furthermore, if £ L°° , then 9 £ L°°([0, T], L°°). 



BOUSSINESQ EQUATIONS 



29 



Proof. We use the notation laid out in Section 2. Let us consider the Galcrkin 
approximation to Pq (or equivalcntly, (5.2)) given by 

(5.3a) (I + a 2 A)d t u n + P n B{u n , u„) = P n P a (6 n e 2 ), 

(5.3b) a t n + P n (V-( Ur A))=O, 

(5.3c) u„(0)=P„u o , 0„(O) = P n o . 

This is a finite dimensional system of ODEs in H n with quadratic polynomial non- 
linearity, and therefore it has a unique local solution in C 1 ([0, T n ), H n ) for some 
T n > 0. Let [0, T*) be the maximal interval of existence and uniqueness of solutions 
to (5.3). We show below that T* = oo for every n. 

Taking the inner product of (1.5c) with 9 n , using Lemma 2.2, and integrating in 
time, we find that for t £ [0, T*), 

(5.4) |e n (t)| = \e n (o)\ < \e \. 

Next, we take the inner product of (1.5a) with u„ and use Lemma 2.1 to find 

Id., ,r 



2 dt 



a 2 l|u ? - 



ie 2 ,u„) < \6 n \\ur, 



(5.5) < \9 W\u n \2+a2\\u n \\2 
Consequently, we have for t £ [0,T*), 

(5.6) \u n (t)\ 2 + « 2 ||u„(t)|| 2 < |u | 2 + a 2 ||u || 2 + t 2 %\ 2 . 

According to (5.4) and (5.6), we see that if T* < oo, then ||u„|| and \8 n \ are both 
bounded in time on [0,T*), and thus the solutions can be continued beyond T*, 
contradicting the definition of T* as the maximal time of existence. Thus, T* = oo 
for all ri £ N. Next, we find bounds on the time derivatives. From now on, we 
work on the interval [0,T], where T was arbitrarily given in the statement of the 
theorem. Using Lemma 2.1, along with (2.16) and (5.4), we have 

< sup |(S(u n ,u n ),P„w)| + sup |(6»„e 2 ,P„w) 



(5.7) 



dt 



< sup \u r , 

l|w||=l 



sup \9 

|w||=l 



w|| — 1 
W 



< |u^|||u ( ")|| + A7 1/2 |6»^| 

< lujllu, 



,\\ + \- 1/2 \8 a \. 

Thanks to (5.6) and (5.7) we obtain that (/ + a 2 A)^^- is uniformly bounded in 
L°°([0, T], V), which implies that ^ is uniformly bounded in L°°([0, T],V), with 



respect to n. Similarly, using Lemma 2.2, we obtain 

dOn 



(5.8) 



dt 



< 



|V2| 



1/2 



which implies that ^jf- is bounded in L°°([0, T},H 2 ) independently of n by virtue 
of (5.4) and (5.6). 

The above bounds allow us to use the Banach-Alaoglu Theorem and the Aubin 
Compactness Theorem (see, e.g., [16, 52]) to extract a subsequence, which we still 
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write as (u n , 9 n ), and elements u and 9, such that 
(5.9a) u„->u strongly in L 2 ([0, T], if), 

(5.9b) u„ u weakly in L 2 ([0, T],V) and weak-* in L°°([0, T], V), 

(5.9c) dun^^ weak-*ini°°([0,T],y), 
at at 

(5.9d) 9 n 9 weakly in L 2 ([0, T] 7 L 2 ) and weak-* in £°°([0, T],L 2 ), 

(5.9c) w cak-* in L°°([Q,T], H' 2 ). 

Next, for arbitrary if and <£>, chosen as in (3.2), let us take the inner product of 
(5.3a) with 4>, and of (5.3b) with <p and integrate in time on [0, T}. After integrating 
by parts several times, we have 
(5.10a) 

- (u„(.s),$'(.s))ds-a 2 / ((u n (s),&(s)))ds + y] / (<u n ,P„^$)ds 

JO JO JO 

= (u„(0),$(0))+a 2 ((u„(0),$(0)))+ / (0„(s)e 2) $(s))ds, 



(5.10b) 

{6 n (s),ip'(s))ds+ [ (6 n u n ,P n V(p)ds = (6 ,<p(0)), 
Jo 

where we have again denoted ' = -4-. We would like to pass to the limit as n — > oo 
to obtain (5.1). The convergence of the linear terms is straight-forward, thanks 
to (5.9). As for the non-linear terms, notice that the convergence in (5.9) is in 
stronger than the convergence in (3.5), and so the convergence of the non- linear 
terms follows just as in the proof of Theorem 3.5 (note that P n <& = $ and P n tp = ip 
for sufficiently large n, due to our choice of test functions in (3.2)). Thus, (u, 9) 
satisfies (5.1). In particular, choosing <f> and $ to have compact support in (0, T), we 
see that the equations of (5.2a) and (5.2b) arc satisfied in the sense of distributions 
in time with values in V' and H~ 2 , respectively. Acting with (5.2a) on <£>, with 
(5.2b) on if, and integrating in time on [to,ii], we find 

(5.11a) 

- [\u{s),&{s))ds-a 2 [\(u{s),&{s)))ds + Y^ I " \u j u,dj$)ds 

J to J to J to 

= (u(t ), *(*o)) + <* 2 ((u(i ), #(*o))) - (u(ti), $(ti)) - a a ((u(*i), $(f x ))) 
+ / \6(s)e 2 ,$(s))ds, 



(5.11b) 

1 (6(s),<p'(s)) ds+ I' (On, V<p) ds = (9(t ), <p(t Q )) - (9(h), ^(h)). 

Temporarily restricting our set of test functions to those which are compactly sup- 
ported in time on [0,T] and considering the case to = and t\ = T, it is easy 
to see from a simple density argument that (u, 9) satisfies the equations of (5.2) 
in the sense of distributions, thanks to (5.11). Next, allowing $(0), and ip(0) to 



BOUSSINESQ EQUATIONS 



31 



be arbitrary, but fixing $(T) = and ip(T) = 0, we act on 4> with (5.2a) and 
on ip with (5.2b) and integrate on [0, T], resulting equations in (5.1). Compar- 
ing (5.1) with (5.11), we find that that 0(0) = 0q in the sense of H^ 1 and that 
(I + a 2 A)u(0) = (I + a 2 A)uq in the sense of V. Inverting I + a 2 A gives u(0) = uo. 
Furthermore, we may send t\ — > to in (5.11b), and use the density of C°°(T 2 ) 
in L 2 (T 2 ), as well as the boundedness of in L°°([0, T], L 2 (T 2 )), to show that 
G C w ([0,T],L 2 (T 2 )). Next, since we have u, G L°°([0,T],V) L 2 ([0, T],V), 
it follows that u G C([0, T], V) by the Sobolev Embedding Theorem. Thus, we have 
shown a weak solution exists in the sense of Definition 5.1. Finally, one may show 
that if O G L°°(T 2 ), then G L°°([0,T],L°°), by following Step 4 of the proof of 
Theorem 3.5 line-by-line. □ 

Theorem 5.4 (Uniqueness for the 2D Voigt model). Let T > be arbitrary. 
Suppose Uo G T>(A) and 9q G L 2 (T 2 ). Then there exists a unique solution to (1.5) 
in the sense of Definition 5.1. Furthermore, it holds that u G L°°([0,T],'D(A)). 

Proof. Here, we only sketch the proof, since the ideas a similar to those given 
above. The existence of solutions to (1.5) has already been established in The- 
orem 5.3. Thanks to the hypothesis uo G T> (A), it is straight-forward to show 
that u G C([0,T],T>(A)) using, e.g., the methods of [34] and similarly that M G 
L 2 ([0,T], H^ 1 ). One can the prove the uniqueness of solutions by following the 
proof of Theorem 3.7 almost line by line. Only some slight modifications to the 
handling of the terms involving ||u||, and in using the parameter a 2 rather than v 
is needed. □ 

Theorem 5.5 (Convergence as a 4 0). Given initial data (uo,0o) £ (H 3 (T 2 ) D 
V) x H 3 (T 2 ), and (ug^) G (H 3 (T 2 ) n V) x H 3 (T 2 ), let (u,6) and (u a ,(9 Q ) be 
the corresponding solutions to the problems Pq and Pq , respectively. Choose an 
arbitrary T G (0,T max ), where T rnax is the maximal time for which a solution to 
the problem P® exists and is unique. Suppose that Uq — > Uo in V and 9q — > 9q in 
L 2 (T 2 ). Then u a -> u in L 2 ([0,T],y) and 9 a -> 9 in L 2 ([0, T], L 2 (T 2 )). 

Proof. Here, for simplicity, we only work formally, but note that the results can 
be made rigorous by using the techniques discussed above. Under the hypotheses 
on the initial conditions, it was proven in [15] that there exists a time T > 
and a unique (u,0) G C([0, T],H 3 (T 2 ) n V) x C([0, T],H 3 (T 2 ) n V) solving the 
problem Pq , (in particular, it holds that T max > 0). Thanks to Theorem (5.3), we 
know that there also exists a unique solution to the problem Pq , namely (u a , 9 a ) G 
C([0, T], U)xC([0, T], L 2 (T 2 )). Subtracting the corresponding equations from these 
problems (that is, (1.1) with v = k = and (1.5)) yields 

(5.12a) 




B(u - u", u) - B(u a 7 u - u Q ) + P a ((9 a - 9)e 2 ), 



((u - u") • V)0 - (u Q • V)(0 - 9 a ). 
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Let us take the inner product of (5.12a) with u a -u and of (5.12b) with 9 a — 8, 
and add the results. After integrating by parts and rearranging the terms, we find 

(5.13) ~ (a 2 ||u- u«f + |u - u«| 2 + 6 a \ 2 ) 

= -(B(u - u", u), u Q - u) + ((9 a - 0)e 2 , u a - u) 
- (((u - u Q ) • V)6, d a -0)- a 2 (Au t , u - u Q ) 

< i|Vuj| L oo|u-u a | 2 + |r-0||u a -u| 
||V»|U~|u - u a \\9 a - 9\ - a 2 (Au t ,u - u a ) 

< if (|u - u Q | 2 + \6 a - 8\ 2 ) - a 2 (Auj, u u Q ) , 

where we have used Young's inequality and the fact that ||Vu||i,«>, ||V0||l°c < oo. 
It remains to estimate the integral on the left-hand side of the equality. Using the 
fact that (u, 0) satisfies (1.1c), we have 

(5.14) -a 2 (Au,,u-u a ) 

= -a 2 (A[-u- Vu- Vp + (9e 2 ],u-u Q ) 

= a 2 (Au • Vu + 2(Vu ■ V)Vu + u • VAu - A9e 2 , u - u Q ) 

< Ca 2 |Au|||Vu|| L oc| u - u Q | + ||u|| L ==|VAu||u - u"| + |A0||u- u Q | 

< Ca 2 ||u|| ff2 ||u|| ff 3lu-u a | + C'||u||^ 2 ||u|| £r 3|u-u Q l + a||0|| if2 |u-u a! l 

< a 2 if|u- u Q |. 

For the second equality, we used (2.18). Combining (5.13) with (5.14) and using 
Gronwall's inequality yields 

a 2 \\u{t) - u a {t)\\ 2 + |u(i) - u a (t)\ 2 + \9{t) - 9 a (t)\ 2 

(5.15) < C (a 2 ||u - <|| 2 + |uo - <| 2 + \6 Q - e«\ 2 ) e *(» 2 +«)'. 

Thus, if Ug — > u in V and 6$ — > 9 in L 2 (T 2 ), as a — > 0, (in particular, if 
u# = u and 6% = 6 for all a > 0), then u" -> u in L°°{[0, T],V) and 9 a -> 9 in 
L°°{[0,T],L 2 ),&sa^0. □ 

Theorem 5.6 (Blow-up criterion). With the same notation and assumptions of 
Theorem 5. 5, suppose that for some < oo , we have 

(5.16) sup limsupa 2 ||u"(<)j| 2 > 0. 

te[o.T») 

Then the solutions to Pq become singular in the time interval [0, T*). 

Proof. To get a contradiction, suppose that (u, 9) stays bounded in H 3 n V(T 2 ) x 
-ff 3 (T 2 ) but that (5.16). Taking the inner product of the momentum equation with 
u a and integrating, we find 

a 2 \\u a (t)\\ 2 + \u a (t)\ 2 = a 2 \\u \\ 2 + |u | 2 + 2 f (9 a (s)e 2 , u a (s)) da. 

Jo 

Taking the limsup as a — > + then by virtue of Theorem 5.3 and Theorem 5.5 we 
have 

(5.17) limsupa 2 ||u a (t)|| 2 + |u(t)| 2 = |u | 2 + 2 f {6(s)e 2 , u(«)) da. 

Q-S.Q+ JO 
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However, given the hypotheses on the initial data, and the well-posedness results 
of [15], it is straight- forward to prove the following energy equality: 

|u(t)| 2 Hu | 2 + 2 [ (6(s),u(s))ds, 
Jo 

so that (5.17) contradicts (5.16). □ 
6. The 3D Boussinesq-Voigt Equations 

We now briefly outline an extension of the previous results to the case of the 
three dimensional Boussinesq-Voigt equations. The details are very similar to the 
2D case, so we only prove formal a priori estimates. In order to control the higher- 
order derivatives, we add a diffusion term to the transport equation. This approach 
is similar to that used in [12, 13, 33, 34] to prove global well-posedness for two 
Voigt-regularizations of the 3D MHD equations. We consider the following system, 
written in functional form, which we refer to as Pq k . 

(6.1a) (I + a 2 A)^+B(u,u) = P CT (0e 3 ), 

at 

(6.1b) 4- + B ( u > 9) = kA6, 

at 

(6.1c) u(0) = u , 9(0) = O . 

Remark 6.1. Note that one could also control the higher-order derivatives by adding 
the Voigt-term ~f3 2 A-^9, /3 > to the left-hand side of (6.1b) to allow for case 
when k == 0. Although the resulting regularized system is well-posed, for the 
sake of brevity, we do not pursue this type of additional Voigt-regularization here. 
However, a similar idea has been investigated in the context of the MHD equations 
in [34] (cf. [33]), and also in [13]. 

Definition 6.2. Let u e V D H 3 (T 3 ), 9 g L 2 (T 3 ). For a given T > 0, we say 
that (u, 9) is a solution to the problem Pq k (in three-dimensions) on the interval 
[0, T] if it satisfies (6.1a) in the sense of L 2 ([0,T],V) and (6.1b) in the sense of 
L 2 ([0,T],iT- 1 ). Furthermore, u g C([0, T], V f) H 3 ), ^ g L°°([0,T],V(A)) H 

L 2 ([o,T],vnH 3 ), 9 g L 2 ([o,niT 1 )nc tB ([o,r] J z 2 ) and ft e L 2 ([0,T],H- 1 ). 

Theorem 6.3. Let u g V D H 3 (T 3 ), 9 g L 2 (T 3 ), and let T > be arbitrary. 
Then there exists a solution to (6.1) in the sense of Definition 6.2. Furthermore, if 
9„ g L'P(T 3 ) for some p g [2,oo], then 9 g L°°([0, T],L P ). In the case 9 g L°°(T 3 ), 
the solution is unique. 

Proof. As mentioned above, we only establish formal a priori estimates here. Sup- 
pose for a moment that 9q g L P (T 3 ). Formally taking the inner product of (6.1b) 
with \9\ p ~ 2 9, p g [2,oo), we find as above that 

(6.2) lS± m i P + K (p-l) f |W| 2 |0r 2 dx = O. 

P at J T 3 

Dropping the term involving k, integrating in time, and sending p — > oo, we find 

for all p g [2, oo], 

\\9(t)\\ LP < \\0oWlp. 
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On the other hand, setting p = 2 in (6.2) and integrating in time, we find 

\9(t)\ 2 + 2 K f \\e(s)\\ 2 ds<\9 \ 2 . 
Jo 

Next, following similar steps as in the derivation of (5.5) and (5.6), we find for 

te[o,T], 

|u(t)| 2 + a 2 ||u(i)|| 2 < |u | 2 + « 2 ||u || 2 + T 2 \9 \ 2 := (K aA ) 2 . 
By formally taking the inner product of (5.2a) with Au, we find 

™(||u|| a + a 2 \Au\ 2 ) = -(B(u, u), An) + (9e 3 , An) 

< C\\u\\\Au\ 2 + \9\\Au\ < CK aA a- l \Au\ 2 + \9 Q \\Au\ 

<C(l+K aiia - 1 )\Au\ 2 + \0 o \ 2 . 

Using Gronwall's inequality, we obtain a constant K a ,2 = (|| u o|| 2 + a 2 |yluo| 2 + 
2 -^-)e CoT , where c := Ca~ 2 + K^a^, such that \\u{t)\\ 2 + a 2 \Au(t)\ 2 < K a , 2 
for a.e. t £ [0,T]. Next, we formally take the inner product of (5.2a) with A 2 u 
(recalling that, in the periodic case, — A = A), to find 

~(\Au\ 2 + a 2 \\Au\\ 2 ) = 5>^u, A 2 u) + (0e 3 , A 2 u) 

3 3 3 

= - ^(Au j ■ dju, Au) -2^2 (diU j ■ d l d j u, Au) - ^(m^/u, Au) 
-((fle 3 ,Au)) 

(6.3) < C\Au\ 2 \\Au\\ + \\9\\\\Au\\ < C {{K aa f + \\Auf + \\9\\ 2 ) . 

These a priori estimates can be used to form a rigorous argument as follows. 
In the case p = 2, existence can be proven by using, e.g., the Galerkin method 
as in the proof of theorem (5.3), substituting the a priori estimates established in 
this section as necessary. Passage to the limit, estimates on the time derivatives, 
and the continuity properties of Definition 6.2 can be established using similar 
ideas to those used in the proofs of some of the previous theorems. For the case 
9q £ _L P ( T 3 ), with p £ (2, oo), we begin by smoothing 9q (e.g., by convolving it with 
a mollificr) to get smooth functions 9^ for each e > 0, which converge to 9q as e —> 
in several relevant norms. Clearly #g £ L 2 (T 3 ). Thus, thanks to the existence of 
solutions for the p = 2 case, there exists a solution (u e , 9 e ) to (6.1) with initial data 
(uo,0§) such that 9 e £ L°°([0,T],L 2 ) and u e £ C([0, T], V C\ H 3 ). One may then 
show that 9 e £ L°°([0, T], H 2 ), e.g., by using the Galerkin method and deriving 
straight-forward a priori estimates on higher derivatives. Since in three dimensions 
L°°([0,T],H 2 ) is an algebra, it follows that \9 t \P- 2 9 e £ L°°([0, T],H 2 ), so that the 
above a priori estimates can be established rigorously for (u e ,9 e ). Furthermore, 
the bounds can be made to be independent of e. Standard arguments show that 
one can extract subsequences of (u e , 9 C ) which converge in several relevant norms as 
e — > to a solution (u, 9) of (6.1) corresponding to initial data (uo, 0q). Taking the 
limit as e — > 0, one may show that 9 £ L°°([0, T], L p ). Finally, in the case p = oo, 
one may employ, e.g., the Hopf-Stampacchia technique used in Step 4 of the proof 
of Theorem 3.5. 
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With the above a priori estimates formally established, we now show the unique- 
ness of solutions to (6.1) under the additional hypothesis that 9o £ L°°(T 3 ). Letting 
(ui, 9i) and (112, 82) be two solutions to (6.1) with initial data (uo, 8q). Let us write 
u := Ui — U2 and 8 = Q\—Qi- As in the proof of Theorem 3.7 we also write £ = A~ 1 8 
and £i = A~ 1 8i, i = 1,2, subject to the side condition J T3 £idx = for i = 1,2. 
Following nearly identical steps to the derivation of (3.13), we find, 

(6-4) ~(\u\ 2 + a 2 \\uf) < C a>K K a!l (\u\ 2 + a 2 \\uf + ||£|| 2 ). 

Similarly, following nearly identical steps to the derivation of (3.12b), we find, 

I|||£|| 2 + K |A£| 2 = (uA6,Ve) +]T(d J u2,V^ 

(6-5) < |u|||^|k~||f|| + l|u 2 |M|£|| 2 < A Q . 4 (|S| 2 + llfll 2 ), 

where K a ^ := max{2|j#o|!L°°, K a ^} < 00. Uniqueness now follows by adding (6.4) 
and (6.5) and using Gronwall's inequality. □ 

7. Appendix 

Wc prove the inequality (2.17). The proof is based on the proof of the Brezis- 
Gallouet inequality [7] and follows almost line- by-line the proof given in [10]. For 
w e T> (A) : let us write 

W = °kW k 

kez 2 \(o,o) 

where are the (normalized) eigenfunctions of A (see Section 2) and := 
(w,Wk). Choose M = (e 1 ^ 174 — l) 1 / 2 for a given e > 0, sufficiently small so 
that M > 1. Wc have 

iiw|| L =o < kki = i° k i + Y, i° k i 

keZ 2 \(0,0) 0<|k|<M |k|>A/ 

\- (1 + lkl 2 ) 1 / 2 , \ ^ (l + |k| 2 ) , 

- 2. (1 + | U | 2) i /2 l a kl+ ^ (iT|kp) |ak| 

0<|k|<M v 1 1 ' |k|>M v 1 1 ' 




y — L - 

^ (l + |k| 2 ) 



1/2 



1/2 

y — ± - 

1/2 / ^ v 1/2 

ax \ 



\J\x\<M (1 + X 2 ) J \J\ X \>M (1 + |x| 2 ) 2 _ 

= C||w||7rlog(l + Af 2 ) + ClAwl^^ 
-cfllwHe-^ + IAwIe- 1 /' 174 
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